#!/bin/sh
cd ${0%/*} || exit 1    # run from this directory


# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions


runApplication blockMesh
runApplication topoSet
runApplication splitMeshRegions -cellZones -overwrite
#runApplication setFields


# remove fluid fields from solid regions (important for post-processing)
for i in wallad2
do
   rm -f 0*/$i/{mut,alphat,epsilon,k,U,p_rgh}
done

for i in topAir wallad2
do
   changeDictionary -region $i > log.changeDictionary.$i 2>&1
done
for i in wallad1
do
   rm -f 0*/$i/{mut,alphat,epsilon,k,U,p_rgh}
done

for i in topAir wallad1
do
   changeDictionary -region $i > log.changeDictionary.$i 2>&1
done

for i in wall
do
   rm -f 0*/$i/{mut,alphat,epsilon,k,U,p_rgh}
done

for i in topAir wall
do
   changeDictionary -region $i > log.changeDictionary.$i 2>&1
done



#-- Run on single processor
runApplication `getApplication`

## Decompose
#runApplication decomposePar -allRegions
#
## Run
#runParallel `getApplication` 4
#
## Reconstruct
#runApplication reconstructPar -allRegions


echo
echo "creating files for paraview post-processing"
echo
paraFoam -touchAll

# ----------------------------------------------------------------- end-of-file
